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^ ■ ABSTRACT 

^ ' Context. Luminous elliptical galaxies generally display a rich star cluster system, whose properties provide strong constraints on the 

O physics of galaxy formation and evolution. Star cluster system studies, however, concentrate on galaxies located in nearby or rich 
galaxy clusters. 

QQ ' Aims. We acquired deep B and I images of NGC 1600, a luminous elliptical in a galaxy group to study its star cluster system. The 

images were obtained with the Optical Imager at the Southern Telescope for Astrophysical Research for an exposure time of 1.66hr 
in each filter. 

. Methods. The sample selection incompleteness was assessed as a function of magnitude and image background level. Source counts 

were measured for different elliptical annuli from the centre of NGC 1600, background subtracted, and fitted with a Gaussian function. 
Colour distributions were derived as a function of galactocentric distance for sources measured successfully in both filters. Typical 
ages and metallicities were estimated based on single stellar population models. 

Results. A clear excess of point sources around NGC 1600 was found in relation to the nearby field. The source counts were consistent 
^ ' with a Gaussian distribution typical of other luminous ellipticals. The luminosity function fits provided an estimate of the density of 

' clusters at the different annuli that could be integrated in solid angle, resulting in an estimated total population of Nac - 2850 star 

clusters. This yielded a specific frequency of 5a> - 1.6. The colour distributions show a hint of bimodality, especially at ^ 20 kpc from 
the centre. Clusters in this region may be associated with a ring or shell perturbation. Finally, the star cluster candidates were cross- 
^ ' correlated to discrete X-ray sources and a coincidence rate of ^ 40% was found. These are likely to be globular clusters harboring 

ly^ ' low-mass X-ray binaries. 

! Conclusions. 

' Key words. Galaxies: individual: NGC 1600 - Galaxies: clusters - Galaxies: elliptical and lenticular, cD Galaxies: formation - 

Galaxies: evolution - 

OO 1. Introduction 
O 

• • Star cluster systems reveal significant information about their host galaxy formation and evolution. By far, luminous ellipticals 
. ^ display the most conspicuous cluster systems and have been the popular targets for studies of star cluster populations (Kundu et al 
k>( , 1999, Harris et al 2006). Studies are also concentrated on nearby galaxies, some in the field, but most often in galaxy clusters, such 
. as Virgo and Fornax (Larsen et al 2001, Santiago & Elson 1996, Kissler-Patig et al 1997, Peng et al 2006a, Strader et al 2006). An 
d ' important result of detailed analyses of extragalactic clusters is that clusters in luminous elliptical galaxies present bimodal colour 
distributions. These are often interpreted as evidence for separate episodes of star and cluster formation, often involving galaxy 
merging (Ashman & Zepf 1992, Cote et al 1998; see recent review by Brodie & Strader 2006). Even though most studies refer to 
extragalactic star cluster systems as globular cluster systems, as in fact globular clusters dominate the cluster counts, the discovery 
of new classes of clusters, such as faint fuzzies (Larsen & Brodie 2000, Brodie & Larsen 2002) and diffuse star clusters (Peng et al 
2006b) suggests that a broader terminology should be adopted. 

Galaxies in lower density environments, such as galaxy groups, have also been studied for the presence of extragalactic star 
clusters (i.e., Kundu & Whitmore 2001, Chies-Santos et al 2006). However, many bright elliptical galaxies in groups still lack an 
investigation of their cluster system. One such case is the luminous elliptical galaxy NGC 1600. 

NGC 1600 is the most luminous elliptical in its loose galaxy group, which also contains two other early-types; NGC 1601 
and NGC 1603. Optical surface photometry, both photographic (Barbon et al 1984) and CCD (Mahabal et al 1995) imaging were 
obtained for the galaxy. Near infra-red surface photometry was completed by Rembold et al (2002). These studies revealed a 
flattened and boxy shape. Sivakoff et al (2004) obtained Chandra X-ray images and identified dozens of point sources superimposed 
on the diffuse emission. They tentatively associated the unresolved sources, some of which are classified as ultra-luminous in X-rays, 
with the population of globular clusters (GCs). 

In this paper, we attempt to detect the star cluster system of NGC 1600 using deep imaging with the Southern Telescope for 
Astrophysical Research (SOAR). Our primary goal is to estimate the total number and specific frequency of clusters in NGC 1600, 
to derive their colour distribution and associate the star cluster sample with the X-ray point sources identified by Sivakoff et al 
(2004). 
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Fig. 1. Close up on the southern side of NGC 1600 using the final SOI/SOAR images. North is upwards and East is to the left of the 
images. Left panel: I band; right panel: B band. 

2. Data and Photometry 

Deep images of NGC 1600 were taken with the SOAR telescope in B and I filters. The SOAR Optical Imager (SOI) consists of 
two E2V CCDs, each one with 2k x 4k pixels, covering a field of view of 5.5x5.5 arcminutes. A 2x2 binning was used, yielding a 
detector scale of 0.154 arcsec/pixel. A total of 6 x 1000s exposures (1.67 hr integration time) were taken in slow read-out mode with 
each filter during two nights. Photometric standards were also observed on the same nights, along with bias and flat-field images. 

The individual exposures were bias-subtracted and flat-fielded using standard tasks from the Image Reduction and Analysis 
Facility (IRAF). The SOI I band images are affected by interference that causes a large-scale fringing pattern. These were corrected 
using the IRAF mscred.rmfringe task, which computes the pattern amplitude minimizing the difference between the target image 
and a fringe pattern image (usually blank sky), and corrects the former The original Multiple Extension Format (MEF) files were 
converted into Flexible Image Transport System (FITS) files, again using the IRAF mscred package. The exposures from the same 
filter were aligned and combined into a final image, which was used for sample selection and photometry. Figure[T]shows a close-up 
of the southern side of NGC 1600 of the final I and B combined SOAR/SOI images. A large number of point sources are seen 
superimposed onto the galaxy's diffuse light. These are the star cluster candidates. The northern side of the image is partly affected 
by a strong diffraction spike from a very bright star Also present in the SOI field (but not in Fig. [1]) are the two other galaxies 
belonging to the same group, NGC 1603 and NGC 1601. 

The IRAF daophot package was used to identify point sources and measure their aperture magnitudes. This was completed 
separately in each of the B and I combined images. Sources with a detection confidence interval 2.5cr above the sky background 
were automatically detected with daofind. An aperture of 4 pixels = 0.61 arcsec was adopted to measure the magnitudes in both 
filters using the task phot from the same package. The magnitudes were measured over the combined image used for source selection 
but after the diffuse light from NGC 1600 was modelled and subtracted-off. Isophote fitting with the IRAF task ellipse was applied 
to model the galaxy intensity distribution. This information was then used to generate a model of NGC 1600 with the task bmodel. 
This model was subtracted from the combined image to produce magnitude measurements that were unaffected by the steep gradient 
caused by NGC 1600 light profile. Since a low threshold above the background was adopted in daofind, a visual inspection of the 
model subtracted image revealed an insignificant amount of extra cluster candidates relative to the authomatic detection. This is 
important, because it is far easier to compute photometric incompleteness in a sample selected by an automatic procedure rather 
than by eye. 

Photometric calibration was based on 36 magnitude measurements in each filter for 7 standard stars taken from Landolt (1992) 
and observed at different airmasses. The magnitudes were measured in the same aperture as the point sources in NGC 1600. The 
calibration fields used for this purpose were SA 95 and PG0231H-051 listed in Landolt (1992). The calibration transformation from 
instrumental magnitudes into Johnson-Cousins system included a linear {B - 1) colour term, a linear airmass term, and a zero-point. 
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Fig. 2. Panel a): Total source counts as a function of calibrated B magnitudes in different elliptical annuli centred on NGC 1600. The 
solid line shows the source counts in the field region, away from NGC 1600 and from the other galaxies in its group. The semi-major 
axis range (in pixels) of the different annuU is given in the lower-right. Panels b,c): The symbols are the field-subtracted counts in 
the same annuli as in panel a. The middle lines are Gaussian fits to the counts, using Eqs. (1) and (2). The lines immediately above 
and below of each fit are Gaussian functions with a 1 cr variation in normalization relative to the fitted one. The final fit, shown by 
a long dashed line in panel c, was carried out for the range 500 < Ripix) < 800 



A least squares fitting of the three coefficients was carried out for each filter separately. Saturated stars and deviant measurements 
were eliminated by applying a 2cr-clipping algorithm iteratively to the magnitude residuals, until convergence. In the B band, a 
total of 1 1 measurements were ehminated in the clipping process; for the remaining 25 measurements, the final fit dispersion was 
ctb = 0.04. In the I-band, 16 measurements were used in the final fit, with a dispersion ctj = 0.06. 

After calibration, magnitudes were corrected for foreground extinction in the direction towards NGC 1600. Schlegel et al (1998) 
quoted Afi = 0.191 and Aj = 0.09. Burstein & Heilles (1982) quoted Ab - 0.08. The average of both Ab determinations was taken. 
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yielding Ab - 0.14. Correspondingly, A/ = 0.065. These values were used to correct the calibrated magnitudes and colours, when 
available. 



2.1. Completeness analysis 

Completeness experiments were performed on the combined images using the daophot.addstar task. Artificial stars spanning a wide 
range of B and I magnitudes were added to the combined images. These were created by fitting a mofFat point spread function 
(with j6 = 1.5) fitted to bright, isolated, and non-saturated point sources found in the images themselves. The fraction of recovered 
artificial stars depended not only on their magnitudes but also on their location, and was lower in regions close to the centre of NGC 
1600, where the background was dominated by the galaxy light. The effect of NGC 1603 and NGC 1601 on sampling completeness 
in the B and I image was also modelled. The background level was locally contaminated by a strong diffraction spike towards the 
north of NGC 1600. This region was masked during the analysis, in addition to the saturated regions in NGC 1600, NGC 1603 (in I 
only), and NGC 1601 (in both B and I images). 

From the completeness experiments, the completeness function could be modelled as a bivariate function of magnitude and 
background level. This was achieved as follows: for each point source found in a given filter, we first determined, by cubic spline 
fits to the completeness vs. magnitude relation, the completeness fraction at its measured magnitude for different sky levels. This 
yielded a completeness versus sky level relation at the source magnitude. Another cubic spline was then applied to this relation 
to evaluate completeness at the background level at the location of the source. Each star was then assigned a weight equal to the 
product l/c(B) X 1 /c(/) where c(B) and c(I) are the completeness fractions calculated for the B and I images, respectively. 



3. Luminosity functions 

Source counts as a function of B and I magnitudes were then derived for several independent elliptical rings around NGC 1600. 
These are shown in Figs. |2] and [3] for the samples selected from the B and I images, respectively. The units are sources kpc"^ 
and the plots are in logarithmic scale. Since a bin size of Amag = 0.5 was used, the numbers were doubled to be normalized to 
unit magnitude. Image surface was converted into physical surface assuming that NGC 1600 is at a distance of 63.0 Mpc (taken 
from NASA Extmgalactic Database, A^/iZ^.This corresponds to a distance modulus of (m - M)q - 34.0 and yields a conversion 
of Ikpc = 3.3arcsec - 2\pix. The radial ranges given in the figure panels correspond to the inner and outer semi-major axis of 
each elliptical annulus used. Panel a shows the total counts, including the background sources located away from NGC 1600 and 
the other galaxies in its group. The remaining panels show the excess source counts around NGC 1600 after the field counts were 
subtracted at each elliptical annulus. 

We assume that all excess counts are star clusters, most of them globular clusters, belonging to NGC 1600 and attempt to fit the 
globular cluster luminosity function (GCLF) at the different radii. A Gaussian GCLF was assumed to have the form: 

^(m) — Poe 2<r'- , (1) 

where we assumed that My = -7.4 and cr - 1.4. This is an intermediate value between that quoted by Harris 1991 and that used 
by Larsen et al (2001) based on HST data. Typical GC colours of (B - V) - 0.85 and (V - 7) = 1.10 were adopted when fitting the 
functions to the B and I counts shown in the figures. These are colours of old (r > 10 Grys) single stellar populations with sub-solar 
(Z ^ 0.008) metallicity, according to Bruzual & Chariot (2003). By fixing the turn-over GCLF luminosity and its width, the fits had 
only the normalization, po, as a free parameter The most accurate Gaussian fits are represented by the middle lines in Figs. |2]and[3] 
They were obtained by varying the normalization po until the quantity 

2 _ ^jv (togO- togOpfa)^ 

reached a minimum. In the expression above, was the number of points to be fitted in each region, 0(m) was the Gaussian 

— 1/2 

function given by eq. (1), <i>ohs is the observed source density. The associated density fluctuation, 5log<^obs ^ A^g/^j , was found by 
propagating the Poisson fluctuation in the source number counts. Nobs, in a given region and magnitude bin. 

The upper and lower lines shown in Figs.|2]and[3]represent ±lcr confidence levels in po for the best-fitting Gaussian functions. 
The variation in po was measured with 300 bootstrap realizations in which the source counts were varied according to a Poisson 
distribution and the fit was repeated. In generating these realizations, a Gaussian distribution of shifts in the assumed distance 
modulus and My, with zero mean and standard deviation o-(„-m) =0.1 and ctmy =0.1 were also applied, to simulate uncertainties 
in these quantities. We emphasize that the adopted typical variation of 0.1 mag in distance modulus is larger than the expected 
uncertainties in the Hubble constant and flow models, and that this variation alters our answers by less than the amplitude of the 
Poisson fluctuations. 

These fits in conjunction with the assumed Gaussian model for the GCLF allow us to estimate the size of the star cluster system 
of NGC 1600. 

The Po values and the expected total number of GCs are listed in Tables [T]and|2]for each annulus. The latter is given by 



^ 2. 2. f* ^22 / 

Ncci^iy^i] = - (^2 - a^) I <I>(m) dm - n — {a2 - a^) po cr yln, 

^ U —CO ^ 



(3) 
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Fig. 3. Panel a): Total source counts as a function of calibrated I magnitudes in different elliptical annuli centred on NGC 1600. The 
solid line shows the source counts in the field region, away from NGC 1600 and from the other galaxies in its group. The semi-major 
axis range (in pixels) of the different annuli is given in the lower-right. Panels b,c): The symbols are field-subtracted counts in the 
same annuli as in panel a. The middle lines are Gaussian fits to the counts, using Eqs. (1) and (2). The Unes immediately above and 
below each fit are Gaussian functions with a 1 cr variation in normalization relative to the fitted one. 



where b/a is the axis ratio of the annulus, oi and 02 are the inner and outer semi-major axes, and (S>{m) is the GCLF in units of 
sources kpc"^, as given in Eq. (1). The axis ratio is obtained from the isophote fitting of NGC 1600; it is b/a = 0.65 with little 
variation with position. Given the range [ai, 02] and the fitted normalization density, we can derive the value of A^cc[«i> «2] with no 
additional free parameters. The uncertainties in po shown in Tables [T] and |2] are the 1 cr spread in this quantity evaluated from the 
results of the 300 bootstrap realizations described earlier The uncertainties in Ncc are propagated from those in po. 

The points marked with an asterisk in Tables [1] and |2] are estimates based on inward extrapolations of the source densities from 
the other bins. They correspond to regions in which the high intensity of the galaxy caused saturation (or a non-linear response 
regime on the CCD) and prevented reliable detection and measurement of magnitudes for a sizable sample. The study of these inner 
regions is however necessary to allow an estimate of the total number of star clusters in NGC 1600. Also shown in the tables is the 
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Table 1. Estimated star cluster surface densities at different elliptical annuli from sample in B filter. 



Range (pix) 


NPTS 


NPTS,, 


poikpc ^) 


Noc (range) 


Ngc (cum.) 


0-100 


* 


* 


1.59 ± 0.40 


253.5 ± 63.8 


253.5 ± 63.8 


100-200 


36 


85.4 


0.82 ± 0.21 


392.2 ± 100.4 


645.7 ± 164.2 


200-350 


121 


267.2 


0.60 ± 0.09 


789.2 ± 118.4 


1434.9 ± 282.6 


350-500 


175 


378.5 


0.38 ± 0.05 


772.5 ± 101.6 


2207.4 ± 384.2 


500-800 


413 


761.2 


0.08 ± 0.02 


497.4 ± 124.4 


2704.8 ± 508.6 



Column 1: range in pixels between inner and outer semi-major axes of annulus used for source counts (1 kpc = 21 pixels); Column 2: total 
number of sources found in the annulus; Column 3: total weighted number of sources in annulus; Column 4: normalization density of Gaussian 
fit to GCLF at given annulus, as explained in text (in units of kpC^); total expected number of clusters in annulus, obtained as explained in 
the text. 



Table 2. Estimated star cluster surface densities at different elliptical annuli from sample in 1 filter 



Range (pix) 


NPTS 


NPTS,, 


pdkpc ^) 


Noc (range) 


Nqc (cum.) 


0-150 


* 


* 


2.20 ± 0.44 


789.2 ± 157.8 


789.2 ± 157.8 


150-250 


79 


156.0 


1.47 ±0.30 


937.5 ± 191.3 


1726.7 ± 349.1 


250-350 


155 


323.4 


0.37 ± 0.06 


354.0 ± 57.4 


2080.7 ± 406.5 


350-500 


191 


353.2 


0.27 ± 0.04 


548.9 ±81.3 


2629.6 ± 487.8 


500-650 


163 


357.6 


0.13 ± 0.02 


357.5 ± 55.0 


2987.1 ± 542.8 



Column 1: range in pixels between inner and outer semi-major axes of annulus used for source counts (1 kpc = 21 pixels); Column 2: total 
number of sources found in the annulus; Column 3: total weighted number of sources in annulus; Column 4: normalization density of Gaussian 
fit to GCLF at given annulus, as explained in text (in units of kpc^^)\ total expected number of clusters in annulus, obtained as explained in 
the text. 



cumulative number of clusters. The system size derived from these counts is almost independent of the band used. The expected total 
number of clusters in NGC 1600 is Ncc - (2.85 + 0.50) 10^. Assuming My - -23.1 to be the absolute magnitude of this galaxy, 
we derive Sn - 1.6 + 0.3 for its specific frequency of clusters. This is relatively low for a luminous elliptical such as NGC 1600. 
Dirsch et al 2003 measured 5 at = 5.1 ± 1.2 for the cD in Fornax cluster, NGC 1399. Rhode & Zepf (2004) found that 5/v = 3.5 ± 0.5 
and S N - 3.6 ± 0.6 for NGC 4406 and NGC 4472, respectively, in Virgo. These galaxies are of similar luminosity but located in 
higher density environments than NGC 1600. It is therefore possible that the difference in 5 at is an environmental effect. However, 
Forbes et al (1997) measured 5 a; = 4.3 + 1.1 for the luminous cD group galaxy NGC 5846. Another possibility is that the star cluster 
system of NGC 1600 extends further beyond the edge of the SOAR image field-of-view, even though the field includes all diffuse 
light from the galaxy. We cannot presently exclude this possibility and this S n estimate should be regarded as a lower limit. 

Figure |4] shows the estimated density of stellar clusters belonging to NGC 1600 as a function of galactocentric distance, ex- 
pressed in terms of semi-major axis of the elliptical annulli. The densities are in units of star clusters kpc"^ and were obtained 
by integration of Eq. (1) over all magnitudes at each distance bin. The points shown as stars or crosses are the result of inward 
extrapolation of the measured densities, as discussed above. Figure|4]also shows, for comparison, the surface brightness profiles of 
NGC 1600 for both filters. These profiles were derived from the isophotal fitting of data for NGC 1600 described earlier We note 
that the cluster system follows the light profile of the field stars out to the projected distances probed (R < 800 pix ^ 40 kpc). This 
is in contrast to that observed for several other luminous ellipticals (Harris 1986), although, as afore mentioned, we cannot exclude 
the existence of an envelope of star clusters outside the field studied. 

4. Colour distributions 

Colour distributions are a useful discriminator between the different star cluster sub-populations, and may reflect different star and 
cluster formation events within a galaxy. Figure |5] shows the {B - I) colour distribution for the point sources that have colour 
information. The units are again in sources kpc"^ and the bin size used was A{B - I) - 0.2. Panel a shows total colour counts of 
sources with B < 21 for the same annuli as in Fig. [3] Panels b and c show the field subtracted source counts. 

A spatial colour gradient is visible in the panels, in the sense that redder star clusters tend to be found closer to the centre of 
the host galaxy. The average star cluster colour varies from (B - /) = 2.2 in the innermost region shown in Fig.|5]to {B - I) - 1.9 
in the 350 < R < 500 pix (17 < R < 25 kpc) annulus. This may reflect true changes in metallicity or age, although redenning 
caused by dust may also contribute. Similar trends were observed by Forbes et al (1997) in the central group elliptical NGC 5846, 
and are also common in luminous cluster ellipticals (Lee & Geisler 1993, Harris et al 1998, Larsen et al 2001). There is a hint of 
bimodality in the distributions associated with NGC 1600 both in the inner regions and beyond R ^ 350 pix (R > 17 kpc). The large 
peak in the inner regions of NGC 1600, with (B - I) ^ 2.2 (see panel b), is typical of an old (t > IQGyrs) and fairly metal-rich 
(Z ^ Zq - 0.019) single stellar population according to the models by Bruzual & Chariot (2003) and assuming extinction to be 
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Fig. 4. The points correspond to the star cluster density profile measured using source counts in the I (solid circles) and B bands 
(open circles). The star and the cross in the inner region are linear extrapolations in the figure of the measured densities in I and B, 
respectively. The curves show the surface brightness profiles in I (dashed line) and B (dotted line), in arbitrary units. 



negUgible. Comparable (B - I) colours were found by Grillmair et al (1999) for NGC 1399 and NGC 1404. The smaller peak at 
(B - I) ^ 1.8 is also present in the region 250 < R < 350 pix (12 < R < 16 kpc), again shown in panel b of Fig.|5] This colour 
requires lower metallicity Z < 0.004 for an old population or significantly younger ages, t ^ 2.5 Gyrs, at solar abundance. 

A clearly larger spread in colours, with two likely peaks, is observed in the intermediate annulus 350 < R < 500 pix (17 < R <25 
kpc), shown in panel c of Fig. |5] The redder peak has approximately the same colour, (B - I) ^ 2.0, as an old star cluster with 
Z ^ 0.008. The bluer peak, however, has (B - I) ^ 1.6. According to Fig. IH this region coincides with the position of a bump 
in the light profile of the galaxy, which is more pronounced in the blue. Probably, this is a ring or shell. Its blue colour is also 
consistent with the smooth profile in the near infra-red bands shown by Rembold et al (2002). In fact, Forbes & Thomson (1992) 
identified an outer shell in NGC 1600 in an R band image. It is possible that the larger spread in colours about this distance may be 
associated with this feature and the bluer cluster peak may represent a younger population of star clusters formed at the onset of the 
perturbation in the host galaxy. 
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Fig. 5. Panel a): Total source counts as a function of calibrated (B-I) colours in different elliptical annuli centred on NGC 1600. 
The solid line shows the source counts in the field region, away from any of the galaxies in NGC1600 group. Panels b,cy. Field 
subtracted counts in the same annuU as in panel a), as indicated. 



5. X-ray counterpars 

The sample of star cluster candidates was cross-referenced with the list of 71 discrete X-ray sources from Sivakoff et al (2004). 
Positional matches were achieved on the basis of the equatorial coordinates. The astrometry of the SOI sources assumed a linear 
relation between CCD and on-sky positions, for the CCD plate scale of 0.154 arcsec/pixel. This assumption was tested with the 
known coordinates of the 3 galaxies imaged, NGC 1600, NGC 1603, and NGC 1601. The differences between their equatorial 
coordinates, as listed in the NASA Extragalactic Database (NED), was compared with the expected differences based on the adopted 
linear transformation from their CCD (x,y) positions. Discrepancies of < 1.5 arcsec were found in all cases, both in decUnation and 
right ascension. 

In attempting to identify the cluster candidates with the X-ray sources we permitted residual offsets in our coordinates relative 
to those hsted by Sivakoff et al (2004). The approach was to start with the null-hypothesis that several optical/X-ray coincidences 
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Fig. 6. On-sky distribution of optical and X-ray sources in NGC1600's group. The small dots are star cluster candidates selected in 
this work. Only sources with magnitudes and colours in the ranges where an excess is seen relative to the field are are plotted. The 
large circles are the discrete X-ray sources from Sivakoff et al (2004). The crosses indicate matched optical/X-ray sources. 



should be found, and then search for the offsets in right ascension and declination, within 3 arcsec from our nominal positions, that 
maximized the number of matches between the two samples. Only optical sources, whose magnitudes and colours were consistent 
with those of star clusters in NGC 1600 (or in its group), were used in this positional matching process. 

Of the 71 discrete X-ray sources identified by Sivakoff et al (2004), a total of 49 were located inside the SOI field limits. 
Sixteen optical/X-ray coincidences were found, and the offsets in equatorial coordinates that resulted in this number of matches 
were 1.75arcsec in right ascension and 0.25arcsec in declination. Figure |6] shows the on-sky distribution of the optical and X- 
ray sources. The two clumps of optical sources correspond to NGC 1600 and NGC 1603. Eleven of the coincidences are clearly 
associated with NGC 1600, being within 40 kpc (132 arcsec) of its centre. A total of 28 X-ray sources are in this same region, 
yielding a coincidence rate of 39%. This may be an underestimate, given the optical sample incompleteness, although there is a 
known trend towards optical/X-ray coincidences being related to bright sources in the optical (Chies-Santos et al 2006, Kundu et al 
2007, Sivakoff et al 2007). For example, Kundu et al (2007) studied 5 early-type galaxies using HST data, and found that the vast 
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Table 3. Sample of Optical/X-ray coincidences. 
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Column 1: source identification number given by Table 2 of Sivakoff et al (2004); Columns 2 and 3: equatorial coordinates (J2000) of X-ray 
sources, given by Sivakoff et al (2004); column 4: X-ray luminosity in units of 10'** ergs/s; columns 5 and 6: calibrated B magnitude and (B-I) 
colour from this work. 



majority of clusters harbouring low-mass X-ray binaries (LMXBs) have M/ < -9. Similarly, most GCs with a LMXB in the sample 
of 11 luminous elliptical and SO galaxies in Virgo studied by Sivakoff et al (2007) have M, < -9. Due to the detection limits of the 
SOAR photometry, we note that all 16 matched sources in this work have Mi < -10 (or My < -8.8) if they are at the distance of 
the NGC 1600 group. 

Kundu et al (2007) and Sivakoff et al (2007) also found that GCs with LMXBs had a strong tendency to belong to the red peak 
of the colour distribution. In the case of NGC 1 399, Kundu et al (2007) found that all clusters with LXMBs had (B-I)> 1 .7. In our 
NGC 1600 sample, we found 4 optical/X-ray matches with bluer {B - 1) values. The remaining 7 objects have redder colours. Since 
the sample was small and we are unable to confirm that individual objects are star clusters in NGC 1600, it is currently impossible 
to assess if this apparent discrepancy is significant. 

Table[3]lists the 16 sources found to coincide with X-ray discrete sources; both optical and X-ray properties, the latter taken from 
Sivakoff et al (2004), are listed. The top 11 entries are those of sources located around NGC 1600. Inspection of Table [3] does not 
reveal any strong systematics between X-ray luminosity and optical colour or magnitude, although the range in the latter quantity is 
small. Only one of the sources listed in the table was previously related to an optical counterpart by Sivakoff et al (2004): this was 
source 40, which is close to the centre of the companion galaxy NGC 1603. All sources matched but one, source 3, are consistent 
with being point sources, according to these authors. Source 3 is close to the centre of NGC 1600. 

6. Conclusions 

Deep B and I photometry of the galaxy group of NGC 1600 has been obtained using the SOAR telescope. Several thousand point 
sources were found in the field for both filters. There was a significant excess of point sources around the luminous elliptical NGC 
1600, which we associated with its star cluster system. Magnitudes and colours were measured and sample incompleteness was also 
quantified. The star cluster luminosity functions at different elliptical annuli around NGC 1600 were all consistent with a Gaussian 
of typical average luminosity and dispersion. Assuming this functional form, the total surface density of star clusters, corrected 
for sampling incompleteness, was obtained as a function of galactocentric distance. The density distribution was similar to the 
underlying light distribution of the galaxy. A total of Ncc = 2850 ± 500 star clusters were estimated to exist within the central 40 
kpc of NGC 1600, yielding a specific frequency of Sn = 1.6 ± 0.3 in this region. Since the cluster system may extend beyond the 
diffuse galaxy light and the limits of the field studied, this estimate should be considered as a lower limit. 

The colour distribution had a hint of bimodality, especially at the centre of NGC 1600 and at an intermediate region ^ 20 kpc 
away. The reddest and most significant peak at the centre is consistent with a t > 10 Gyr old metal-rich (Z ^ Z©) star cluster 
population. The bluer peak, which also dominates the cluster counts at 10-15 kpc, is either an old population with 1/5 of the 
solar abundance or a 2.5 Gyr old population of solar metallicity. 

The colour distribution became particularly wide and bimodal at ^ 20 kpc from the centre of NGC 1600. This region coincided 
with a bump in the surface brightness profile of the host galaxy, which is clear in the B filter image and completely absent in the 
near infra-red. This is probably a ring or shell, as previously suggested by Forbes & Thomson (1992). This perturbation may have 
been accompanied by star and cluster formation. The bluer colours found in the cluster sample in this region may be related to such 
a star formation event, which would indicate a relatively young star cluster population. 

Finally, the optical cluster sample was cross-correlated with the discrete X-ray sources from Sivakoff et al (2004). The coin- 
cidence rate within 40 kpc of the galaxy centre, where an excess of point sources in the optical is clearly seen, is ^ 40%. The 
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optical/X-ray coincidences vary widely in X-ray luminosity and optical colours, but due to the magnitude limit of the SOAR im- 
age they all correspond to optically luminous star clusters (Mj < -10). Since previous studies, based on Chandra and HST data, 
indicated that clusters with LMXBs are preferentially luminous, it is possible that the sample of optical/X-ray matches that we have 
found is a reasonably complete census of GCs harbouring X-ray sources in NGC 1600. Except for one source, all matched sources 

are consistent with being unresolved in the Chandra images. 
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